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Abstract 

VT— symmetrization of quantum graphs is proposed as an innovation where 
an adjustable, tunable nonlocality is admitted. The proposal generalizes 
the VT— symmetric square-well models of Ref. [1] (with real spectrum and 
with a variable fundamental length 9) which are reclassified as the most 
elementary quantum q— pointed-star graphs with minimal q = 2. Their 
equilateral q — 3,4, . . . generalizations are considered, with interactions at- 
tached to the vertices. Runge-Kutta discretization of coordinates simplifies 
the quantitative analysis by reducing our graphs to star-shaped lattices of 
N = qK + 1 points. The resulting bound-state spectra are found real in an 
iV— independent interval of couplings A G (—1, 1). Inside this interval the set 
of closed-form metrics 0^ (A) is constructed, defining independent eligible 
local (at j = 0) or increasingly nonlocal (at j = 1,2,.. .) inner products in 
the respective physical Hilbert spaces of states Hj(X). In this way each 
graph is assigned a menu of non-equivalent, optional probabilistic quantum 
interpretations. 



1 Introduction 



Many nontrivial quantum systems are described via a simplified effective 
model. Vibrational excitations of fields, nuclei or molecules may often be 
represented, for example, by artificial models where a single real or virtual 
(quasi)particle moves along a suitable one-dimensional trajectory, finite or 
infinite. In paper I [I] we even analyzed a family of models where this trajec- 
tory has further been replaced, in the so called Runge-Kutta approximation, 
by a finite lattice of points. 

Specific difficulties may survive even after a drastic reduction of the num- 
ber of degrees of freedom. Typically, a quasi-one-dimensional narrow-tube 
trajectory may happen to be curved (causing the emergence of bound states 
[2]) or twisted (returning these bound states to the free-motion continuum 
again [3]). Other pathologies may emerge when the (quasi)particle moves 
along a topologically nontrivial waveguide. In the simplest thin-tube real- 
ization of the latter scenario one speaks, in general, about the motion of a 
quantum (quasi)particle along a graph, i.e., along a system of one- dimensional 
free-motion trajectories (called "edges" of the graph) connected at the so 
called "vertices" of the graph (various nontrivial interactions could be ad- 
mitted at these points). 

A purely phenomenological motivation of interest in quantum graphs has 
originally emerged in quantum chemistry where the edges were identified 
with the bonds between atoms in a larger organic molecule along which the 
electrons might move almost freely [I]. Soon, a more abstract appeal of 
quantum graphs prevailed offering a nontrivial quantitative picture of quan- 
tum dynamics in many arrangements ranging from the Y-shaped tree up to 
fractal trajectories [5]. The recent proceedings [6] can be cited as a source 
of updated information about the current state of the art. More than 700 
pages of predominantly mathematically oriented reviews still incorporate a 
few physics-centered summaries of potentially appealing phenomenological 
consequences and applications of the theory. Pars pro toto we could point 
out refs. [H |S]) putting more emphasis on physics and listing many related 
references. Today, the use of quantum graphs ranges from the analysis of 
photonic crystals up to the studies of thin wires and waveguides and other 
mesoscopic devices produced by sophisticated nanotechnologies. On theo- 
retical level quantum graphs are increasingly popular as formal structures 
testing field theory [9] or describing certain important phenomena in solid- 
state physics [10]. Multiple concrete models serve as a laboratory of our 
understanding of systems with constraints [UJ. Last but not least one finds 
quantum graphs used as benchmark systems in quantum chaos [12] and/or 
random walks [15] . 

The incessant transfer of the quantum-graph idea from its original, purely 
descriptive role to a more abstract theoretical framework may be expected to 
continue. An illustration of the emergence of new tendencies in this field may 
be seen, e.g., in the complexified, non-Hermitian boundary-supported inter- 
actions as studied, in the context of fully realistic three-dimensional lattice 
models, in Ref. [14J. These tendencies grew from origins which may be traced 
back to a few papers by Hatano and Nelson [J5] and by Feinberg and Zee [H] 
as well as to a number of more recent studies rooted not only in solid-state 
physics [10J but also, say, in nuclear physics [17] . field theory [IB] or cosmol- 
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ogy [19]. An intensification of interest in all of these non-Hermitian quantum 
systems with real spectra occurred, in particular, after the publication of 
influential letter [20] where the very special form of non-Hermiticity called 
PT— symmetry of Hamiltonians (which is to be explained below) has been 
promoted as an unexpectedly productive heuristic principle. It found many 
concrete applications reported, e.g., in proceedings of several dedicated con- 
ferences [2T] . Virtually all of these studies may be characterized as a search 
for a new point of optimal balance between the mathematical requirements of 
simplicity (and, in particular, of constructive tractability of physical models) 
and the natural requirements of dynamical and phenomenological relevance 
of new models. For illustration let us mention just the recent proposals of the 
tests of quantum brachistochrones (where the recent letter [22] summarizes 
the existing theoretical proposals) or of measurements over certain anoma- 
lous scattering systems [23]. The quickly developing discussion of possible 
measurable effects involving VT— symmetric systems in quantum optics [21] 
(and, perhaps, in quantum gravity etc [22]) must also certainly be mentioned 
here. 

The emergence of all of these new theoretical ideas motivated also our 
present work. Their multisided applicability persuaded us not only about an 
undeniable phenomenological appeal and relevance of non-Hermitian interac- 
tion models but also about the promising tractability and feasibility of many 
of their computational and constructive aspects. In what follows we shall 
propose and study, therefore, a schematic though still nontrivial quantum- 
graph models based on a non-Hermitian form of interaction supported, as 
usual, just by certain vertices of the given graph. 

The text will start from a concise outline of the inspiration and origins 
incorporating simple square-well models reviewed in Section [21 and their ele- 
mentary Hermitian star-shaped discrete quantum graph generalizations pro- 
posed in Section [3J The key ideas of our present innovations will be then 
listed in Section H] followed by Section [5] where the main necessary property 
of our quantum graph models, viz., the reality of their spectra will be demon- 
strated. Section [6] will then be devoted to the presentation and explanation 
of the core of our message, viz., to the description of a few first nontrivial ex- 
amples of nonlocal VT — symmetric quantum graphs. Finally, our concluding 
remarks will be collected in Section [7J 



2 Square-well Schrodinger equations 
2.1 Runge-Kutta discretization 

Let us start our considerations from the most common ordinary differential 
Schrodinger equation for bound states in a square well, 

- |U(0 = ^(0, V(±£)=o (i) 

and review a few results obtained for various modifications and perturbations 
of this model in the recent literature. First, let us mention the study [26] 
where the addition of a "sufficiently small" potential V(£) has been shown 
to leave the spectrum real (i.e., in principle, observable) irrespectively of the 
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detailed form of function V(£) which can even be allowed complex. This 
result can be perceived as one of the most persuasive rigorous mathemat- 
ical confirmations of the Bender's and Boettcher's conjecture [20] that in 
applied Quantum Mechanics there exists a broad class of complex potentials 
supporting real spectra of bound-state energies. 

Bender et al [27] noticed and emphasized that many Hamiltonians H = 
p 2 + V(£) 7^ H can be characterized by their VT— symmetry, i.e., by the 
property HVT = VTH with parity V and with time reversal T mimicked by 
Hermitian conjugation [28J. Bound states in a few solvable VT— symmetric 
piecewise-constant potentials were studied in Refs. [29]. The correct proba- 
bilistic interpretation of some of these potentials found its first constructive 
formulation in Ref. [3D]. In parallel, an efficient simplification of the under- 
lying mathematics via Runge-Kutta (RK) discretization of coordinates has 
been proposed in Refs. [31] . It was based on the replacement of the interval 
of £ G [-L, L) by its discrete version 



£*: = kh , k = 0, ±1, 
i.e., by the discrete lattice of points 

€-k+i ■ 



±K , h > 



£- 



K 



£_ 2 £_! £ £i£ 



IK 



(2) 



(3) 



In this perspective one has to replace differential Eq. ([T]) by its discrete 
analogue or approximation 



h 2 



(4) 



A clear insight in the formal structure of the square-well eigenvalue problem is 
achieved. With —K < k < K and ^(£±(#+1)) = our difference Schrodinger 
Eq. (jU) reads 
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i.e., it acquires the transparent matrix-diagonalization form. 
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(5) 



2.2 Equivalence to a linear discrete quantum graph 

Let us renumber the linear array (j3J) of AT = 2K + 1 points in a slightly 
unusual manner which emphasizes its left-right symmetry, 



X2K-1 X 2 K-3 ■■■ X 3 Xi X X 2 X A 



X2K-2 X 2 K 



(6) 
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Our Schrodinger Eq. (jSJ) becomes rearranged, 
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ip(x ) 
ijj(x 2 ) 



(7) 

It may be perceived as describing a system which lives on the linear (one 
could also call it V-shaped) graph which consists of two wedges connected in 
the origin. 

The iV— dimensional Hamiltonian with iV = 2K + 1 as it appears in 
Eq. has a block-tridiagonal partitioned matrix structure 
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(8) 



with v = (—1, —1) and = (0,0) being two-dimensional row vectors while 
2 is a number. The rest of the matrix is composed of two-dimensional 



u 



unit matrices I and null-matrices 0. In the light of what has been written 
in Introduction the naive, discrete quantum square-well problem may be 
reinterpreted as one of the simplest quantum graphs, therefore. 



3 Star-shaped discrete quantum graphs 

The example of preceding section may be complemented by a series of its 
generalizations living on q— pointed star graphs with q = 3, 4, . . .. In this new 
context the trivial example (151) + (171) indicates how this generalization can 
be "translated" back into the language of difference or matrix Schrodinger 
equations. Let us now complement this idea by a few concrete examples of 
its implementation. 

3.1 Y-shaped model: q = 3 

The simplest nontrivial discrete realization of a graph with q = 3 may be visu- 
alized as an Y-shaped (or, if you wish, T-shaped) N— point lattice composed 
of three equally long branches whose individual points will be numbered as 
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follows, 



Xn-2 



X 5 x 2 



x 



X 3 



Xn-4 %N-l 



X4 



(9) 



XN- 



XN-3 



This lattice connects the three (framed) endpoints with the central (doubly 
framed) junction at Xq. The simplest version of a quantum system living on 
this graph may/will employ again the RK discretization of the kinetic energy 
(i.e., of the second derivative operator, cf. Eq. pi)). The only exception is 
encountered at xq where our choice of an acceptable matching is more flexible 
(see Ref. [32J ) . For the sake of simplicity we shall postulate 



ip(xi) + ip{x 2 ) j V j){x q ) -wip{xp) 

h 2 



Ei/j(x ) 



(10) 



with a free parameter u = u(q) set equal, say, to 3 at q = 3. In the bound- 
state arrangement this matching condition in the origin must be comple- 
mented by the three "asymptotic" Dirichlet boundary conditions imposed at 
the remote ends of the edges. The bound-state energies will then coincide 
with the eigenvalues of the real and symmetric (3K+ 1)— dimensional matrix 
Hamiltonian with partitioned structure shown in Eq. (jSJ). Wave functions will 
be specified by Schrodinger equation 
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(11) 

From the symmetry (i.e., Hermiticity) of the Hamiltonian one deduces that 
at any integer K the spectrum is real though not necessarily nondegenerate. 
At iV = 4, for example, we get = 2 while E^l = 5/2 =f \/T3/2. 
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3.2 X-shaped model and its star-shaped descendants 
with q > 4 

At q = 4 the lattice-points should be numbered in the same manner as above, 



XN-2 



x N . 



XN-3 



XN-7 



Xq X 2 



X 3 



Xq 



X% 



X.i 



■ ■ ■ X N - 5 



XN-1 



(12) 



XN- 



xn-a 



The extension of this pattern to any positive integer q is obvious. In the 
corresponding Hamiltonian ([8]) we may keep the RK-discretization-related 
scalar parameter u = u(q) variable or equal to its "maximum" u(q) = q 
tractable as natural after embedding of our graph into a sufficiently high- 
dimensional space. We may add that at any u the degeneracy of the spectrum 
will grow with q. For illustration we may use the model with the smallest 
dimensions N = N(q) = q + 1 where the energy eigenvalue E = 2 proves 
(q — 1)— times degenerate. This is easily seen from Eq. (fTTl) and/or from its 
q > 3 generalizations once we put there, tentatively, i/j(xq) = 0. The whole 
set of equations then degenerates to the single constraint Y^j=i ^(^j) — u — 2 
with q — 1 linearly independent eigenvector solutions. 

At the two remaining unknown energies E = -Ei^+i ^ 2 we may normalize 
ip{xo) = 1 and eliminate ip{xj) = 1/(2 — E) at all j > 0. We arrive at the 
elementary Bethe-ansatz-type quadratic secular equation q/(2 — E) = u — E 
giving the two missing roots in closed form, 



2 E hq+1 = 2 + «Tv / (2-«) 2 + 4g, 



(13) 



This is the first nontrivial q— star-graph-spectrum formula which is, of course, 
compatible with its above-mentioned special case computed at u(q) — 3 for 
q = 3. 



4 Innovation: Two changes of perspective 

The message delivered by the examples presented in preceding sections can 
be summarized as a recommendation that the current discrete square-well 
eigenvalue problem with q = 2 can easily be generalized to its q— pointed- 
star analogues with any integer q > 2. Formally these models may be 
characterized by the iV— dimensional partitioned Hamiltonian matrices H^ N ' 
of Eq. ([H]) where we set iV = qK + 1 and use q— dimensional row vectors 
v = (—1, —1, . . . , —1) and = (0, 0, . . . , 0) and q— dimensional unit matri- 
ces and null-matrices / and 0, respectively. Of course, nothing really new 
emerges in such an elementary constructive project which requires just a 
routine application of the well known principles of quantum mechanics. 
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The situation becomes much more exciting when the purely kinetic nature 
of the Hamiltonian of a quantum graph is enriched by an interaction added, 
preferably, at the vertices. For us, this option opened a way toward two 
generalizations which will be described in what follows. In essence, they will 
be based on the thorough change 

• of the naively Hermitian nature of the interactions (we shall advocate 
here the transition from the usual real and symmetric interaction ma- 
trices if ( m *) to their asymmetric alternatives preserving the reality of 
the spectrum, cf. paragraph 14. II below) . 

• of the naively realistic assumption of the strict locality of the mod- 
els (this will represent a further development of the idea proposed in 
Ref . p] and briefly recalled in paragraph 14.21 below) . 

4.1 VT— symmetric interactions at vertices 

One of the purposes of our present text is to enrich the picture of dynamics 
of bound states living on quantum graphs via an introduction of certain 
nontrivial interactions at their vertices. In a broader physical context this is 
the project inspired not only by Ref. pQ (on bound states) but also by some of 
our other papers (dealing with scattering). In the language of mathematics, 
the formal connections between these two physical scenarios are quite close, 
especially in the RK discretized models. Thus, although there is no space 
here for a deeper study of the scattering on the PT— symmetric graphs, we 
find it meaningful to mention, briefly, at least some of the possible parallels. 

4.1.1 A brief detour to scattering models 

In our few recent papers on scattering [33J E31 ESJ EE] the introduction of cer- 
tain elementary nearest-neighbor VT — symmetric interactions between RK 
lattice points proved fruitful as a very useful and productive model-building 
principle. Unfortunately, there exist several obstacles for making the analogy 
between the bound- and scattering-state one-dimensional RK-based models 
sufficiently close. Firstly, one must keep in mind that in the scattering sce- 
nario the number N of the RK lattice points must be kept very large or 
infinite. Secondly, the very essence of the arrangement of the scattering 
experiments requires that the interactions themselves should preferably be 
localized very close to the origin [331 EZ]- In contrast, the bound-state ar- 
rangement of Schrodinger equations seems to prefer the transfer of the sup- 
port of interactions to the remote ends of the interval of coordinates. In such 
a case, perceivable technical simplifications were reported not only in the 
one-dimensional continuous-coordinate square- well models (cf. Refs. [38J) 
but also in the realistic three-dimensional discrete-lattice calculations (cf. 
Ref. [EJ). 

This being said, a note on some lattice-based models of scattering may 
still prove approprate. Firstly we could classify them more easily in our 
present graph-based language. The presence of a nearest-neighbor coupling 
will be indicated by the insertion of symbol ^> between the corresponding 
two lattice points. In the scattering-inspired arrangement these points are 
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usually chosen as lying not too far from the origin. In the first step the 
following modification is obtained of the discrete graph of Eq. (jSJ), 



X2K-1 X2K-3 



%3 X\ II XO II X 2 X 4 ••• X 2 K-2 [x~2K 



(14) 



For illustration of a quantum system living on this graph let us recall the 
interaction matrix of Ref. 1 33 1 , 



9 
-9 
... 







(15) 



This real and antisymmetric (i.e., VT— symmetric [33j ) matrix has to be 
added to the purely kinematic discrete square- well Hamiltonian (|S]). Unfor- 
tunately, some of the predictions of this oversimplified model are unphysical 
|39j . In subsequent Ref. [31] another version of VT— symmetric interaction 
has been proposed, therefore. It employed the fully symmetrized localization 
of the nearest-neighbor interactions in the RK graph, 



X2K-1 X 2 K-3 . . . X 3 Xi X <} X 2 X 4 . . . X 2 K-2 X 2 K 



(16) 



leading to the amended interaction matrix 



9 9 .. 

-9 ... 

-9 ... 

... '•• 



(17) 



In Refs. [M] and [36] we further shifted the diamonds ^> (representing the 
localization of interactions) by one step in the lattice and arrived at the next 
graph 



X2L-1 



X2L-3 



X 3 Xi 
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x 2 X A 
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X2L 



yielding the next eligible interaction matrix 
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(19) 



A general pattern emerges clearly. The whole class of interactions can be 
realized via four nonvanishing matrix elements which are not necessarily 
located just in the closest vicinity of the origin. This quadruplet of off- 
diagonal matrix elements is allowed to move away from the origin forming a 
series of descendants of Eq. ( Fl9l) . 
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The main benefit of this series of models defined on RK lattices is three- 
fold. Firstly, their study opens the way toward the unitary scattering systems 
described by the sufficiently elementary VT — symmetric Hamiltonians [33J. 
Secondly, the physical predictions (i.e., the reflection and transmission co- 
efficients) retain the form of closed formulas [34|. Thirdly, these models of 
scattering may find generalizations living on some suitable classes of nontriv- 
ial quantum graphs in the nearest future. 



4.1.2 VT — symmetric bound-state models with q = 2 



In contrast to the scattering scenario where, typically, the matrix in Eq. (I19p 
is infinite-dimensional, the RK version of the bound-state problem may al- 
ways be considered finite-dimensional. Then, the repeatedly shifted symbol 
<0 of the interaction must ultimately reach the ends of the V-shaped graph, 



X2L-1 



%2L- 



x 3 Xi 



X 2 X 4 ... X 2L -2 



%2L 



(20) 



The related exceptional Hamiltonian matrix represents the modified square 
well with a nontrivial VT— symmetric interaction which is localized solely in 
the closest vicinity of the external vertices. The related quantum Hamilto- 
nian acquires the partitioned (K + 1)— dimensional tridiagonal form 
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i.e., at q = 2, 
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The rightmost lowest corner carries all the dependence of the Hamiltonian 
on the coupling (note that we changed its symbol from g to A) . This parallels 
the preferences recommended in Refs. [38J or 
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4.2 Introduction of nonlocality via inner products 

Whenever we declare a matrix [e.g., our Hamiltonian ( 12TI) considered in the 
RK coordinate representation] manifestly non- Hermit ian, we almost always 
have in mind just the non- Hermit icity in the current £2 —representation of 
the Hilbert space. This space may be denoted by the symbol 7v- F ' where the 
superscript stands for the "first" or "friendly" space (cf. also [10]). In this 
space the usual formula 

JV-l 

52riMM*k) :=Mih) (22) 

fe=0 

defines the inner product between any pair of its elements (i.e., finite- or 
infinite-dimensional vectors) ip\ and ip2- In this setting the authors of Ref. jTTJ 
noticed and emphasized that the same Hamiltonian may appear to be Her- 
mitian in another Hilbert spaces Ti.( s > where our choice of the superscript 
stands for the "second" or "subtle" space and where the same set of vec- 
tors is merely assigned the following different, non— £2 inner product using a 
suitable nontrivial "metric" ^ I, 

N-l N-l 

J2J2^ e ^ x ^M = <^i|©l^2> := ((^2) • (23) 

One can always make use of this flexibility of basic definitions, keeping only 
in mind that the standard probabilistic interpretation can solely be assigned 
to a Hamiltonian which is Hermitian (in whatever Hilbert space). In this 
sense the models described by asymmetric real Hamiltonian matrices with 
real spectra do not leave the territory of the standard formalism of quantum 
theory. 

Although the latter idea has thoroughly been explained by several authors 
[T71 HU H21 H3] , some of its key aspects and consequences may be summa- 
rized in two brief sentences. Firstly, we must assume that the spectrum of 
our Hamiltonians H = H(X) remains real in some non-empty interval of the 
measures of their asymmetry A. In the second step we have to introduce an 
invertible operator f2 which maps our £ 2 Hilbert space "HS F ' onto another, 
unitarily non-equivalent "physical" £2 Hilbert space TiS p > which is expected 
unitarily equivalent to the "subtle" physical space endowed with non- 
trivial metric and product (|23|) . 

More thoroughly, both these steps will be explained in section |4.2. II Now 
let us only add that their practical appeal has been well illustrated in nuclear 
physics where 1HS P ^ represented the textbook Hilbert space of nucleons (i.e., 
fermions) while both the auxiliary Hilbert spaces H.( F ' S ^ were identified with 
the spaces of certain artificial, effective "interacting bosons" (cf. Ref. [17] for 
more details). 

In Ref. [T] we also worked with the triplet of spaces 

H (F,S,P) and 

em- 
phasized there the deep technical nontriviality of the construction of the 
necessary metric operator = Q(H) in terms of which our asymmetric, 
non-Hermitian real-matrix representations of the Hamiltonians were made 
Hermitian with respect to the ad hoc inner product (|23|) . In particular, as 
long as we worked in coordinate representation, we made distinction between 
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the models which were local (i.e., where = Qq(H) ^ I remained represented 
by a diagonal matrix) and nonlocal, i.e., characterized by the non-diagonal 
metrics Qi(H) , 6 2 (H) , . . .. Moreover, the most unexpected property of the 
oversimplified models as studied in Ref. [I] has been revealed in the fact that, 
via a suitable renumbering, one could achieve that the j—th metric Qj(H) 
was represented by a very special (2j + 1)— diagonal matrix. 

In our present paper we intend to demonstrate that these results may be 
extended to a broad family of quantum graphs. 



4.2.1 A return to (hidden) Hermiticity of observables 

In the formalism described in Ref. [17] the simple but non-Hermitian Hamil- 
tonian matrix H ^ W defined in Ti^ has been put in correspondence with 
its idealized isospectral partner f) = ty N >(X). The latter operator is defined 
in H {p) and it may be assumed complicated. The correspondence is realized 
by the Dyson map, 

: H -> f) = VtHVt- 1 (24) 

which is, by definition, non- unitary, Q = flW(A) + (l/ft) f . Thus, we are 
allowed to require the Hermiticity of the isospectral partner Hamiltonian, 

f) (JV) (A) = n(A) H {N) (X) n-\X) = [fj (JV) ] f (A) . (25) 

The latter relation can be re-read as a constraint imposed upon the simpler 
operator H = H^ N \X) itself, 

qh^q- 1 = [qh^q- 1 ] 1 = [ft- 1 ] 1 " [#W] + nt. 

In the re-arranged and abbreviated form this relation coincides with the 
condition of a hidden Hermiticity or "quasi-Hermiticity" [TTJHl] of H^ N \\), 

[H^ = &H^e-\ = fttft>O. (26) 

The closest correspondence between metric G of Eq. (123]) and the Dyson map 
Q is established in this manner. 



4.2.2 The reconstruction of the ad hoc metric = 0(H) 

In any VT— symmetric quantum model, i.e., for Hamiltonians with the prop- 
erty W = V HV^ 1 (and with the real spectrum) the correct physical prob- 
abilistic interpretation of bound states must be based on the reconstruction 
of the metric in TiS s \ The matrix elements of this metric may be made 
available as a solution of the linear algebraic system of Eqs. 



N 

\{H*) jk e kn -e jk H kn ]=0, j,n = 1,2,..., N. (27) 



fc=i 



Needless to repeat that the resulting metrics are Hamiltonian-dependent and 

by far not unique in general, = Qj(H), j — 0, 1, Each of them defines 

a new, independent Hermitian conjugation and, hence, a respective inde- 
pendent N— dimensional physical Hilbert space 

H (S) = n (N)_ The knowl . 

edge of the metric is substantial. The parallel availability of the factor 
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Q and of its conjugate flt remains less essential (though note their role 
in section 15.2.21 below). For this reason we shall pay our main attention 
here just to the constructive assignment of one or several alternative metrics 
= Qj (A), j = 0, 1, ... to a given, "prescribed" quantum-graph Hamilto- 
nian H = H^ N \\). 

This project consists of fulfilling two separate subtasks. Firstly, we shall 
search for the metric in the form of a superposition 

Q = ^L.=^ N) +pM N) + --- (28) 

of some suitable Hermitian, sufficiently simple though not necessarily positive 
definite auxiliary components. Secondly, due attention must be paid to the 
positive definiteness of the metric ( 1281) controlled by the appropriate choice of 
parameters (3j. In addition, all of the "pseudometrics" will individually 
be assumed compatible with the Hermiticity condition ( 1261) . 

N 

0, j,n = l,2,...,N, 



k=l 



V = 0, 1, . . . . (29) 

In this manner our ansatz (128]) will specify metrics as superpositions of 
pseudometrics V = P^' = v[ N \V2 N \ ■ ■ ■ which will be required to pos- 
sess a sparse-matrix structure. At q = 2 this idea has been shown pro- 
ductive in Ref. [T|. In our present paper we just extend this recipe to the 
VT— symmetric quantum graphs with q > 3. 



4.2.3 Nonlocal metrics: their sample construction at N = 4 

For quantum systems living on the smallest Y-shaped discrete graph 



x 2 


x 


x 3 




Xi 





(30) 



no space is left for the end-point additional interactions since we do not wish 
that the matching point xq gets involved. Thus, our present N = 4 Y-shaped 
quantum graph will remain purely kinematic. Its spectrum of energies will 
coincide with the eigenvalues of the four-dimensional matrix Hamiltonian 

3 -1 -1 -1 " 
-12 

-1 2 ' (31) 
-10 2 

The latter particular matrix is real, symmetric (i.e., Hermitian) and positive 
definite. These properties (plus its natural commutativity with itself) make 
this matrix eligible as an admissible metric. Further metrics compatible with 
their implicit algebraic definition (|26|) can be sought as arbitrary polynomial 
functions of Hamiltonian (l3~Tj) . 

© (4) = ©{cL...} = co I + ci # (4) (0) + c 2 [tf (4) (0)] 2 + ■ • • . (32) 
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The recipe is quick since the necessary explicit construction of the integer 
powers of the Hamiltonian is straightforward yielding 



12 -5 -5 -5 " 
-5511 

(33) 

-5151 
-5 1 1 5 

etc. Unfortunately, the construction of metrics via Eq. (|32l) cannot be trans- 
ferred to non-Hermitian matrices H ^ W . Another unpleasant feature of the 
metrics sampled by Eq. (133]) lies in their non-sparse, full-matrix form. For 
both of these reasons a return is recommended to the methods of paragraph 
14.2.21 Their results are universal - for example, metric Eq. (133]) appears as 
a special case of formula (138]) [cf. Section [67T1 below] at a = 12, b = —5 and 
f = i = k = \. 



ei 4 ok..., = [tf ,4, (o)] 2 = 



5 The proofs of the reality of energies 

At any integer number q the energy spectrum of our quantum graphs is par- 
tially degenerate at A = 0. This leaves the specification of a complete basis 
ambiguous. Another ambiguity emerges via the non-Dirac metrics 0j 7^ /, 
j = 0, 1, We may construct several alternative, nonequivalent repre- 
sentations of the respective "correct" or "selected" Hilbert space of states 
H W (X) = nf\\), j = 0,1,.... In the respective inner products (123]) 

one encounters mutually nonequivalent metrics 0^(A) sampled in para- 
graph U]23] above. The discovery of such a new freedom of making the choice 
between alternative inner products can be perceived as belonging to the most 
important recent achievements in quantum physics, with impact ranging from 
the new flexibility of the interacting boson models in nuclear physics [T7j and 
from formulations of several new theoretical ideas in quantum mechanics (12] 
up to the emergence of the new classes of phenomenological Lagrangians in 
quantum field theory [UJ where, e.g., the presence of ghosts can successfully 
be eliminated in some cases [45] and where even the concept of integrability 
acquired an updated meaning |43j . The use of the varying non-Dirac metrics 
8^ / also opened the way toward new challenges connected, e.g., with the 
description of bound states in time-dependent systems [26] or in the relativis- 
tic kinematical regime |47j. In some phenomenological models of scattering 
the variability of 6 has been suggested as a guarantee of the causality and/or 
unitarity of the process [33] EH EH EH] • 

In our present treatment of the ambiguity of Q = Qj, j = 0,1, ... we shall 
be guided by the approach of paper \J\. We considered there the standard 
coordinate representation (a: 1 of the metric operator and required that 
a suitable measure of its "nonlocality" [i.e., of its deviation from the Dirac's 
"local" delta function 5(x — x')] should be identified with the postulate of 
the existence of fundamental length 9 which characterizes the physical sys- 
tem in question. The same philosophy will also be accepted in our present 
text. We shall assume that the appeal of the concept of fundamental length 
survives the transition to the discrete-graph Hamiltonians of any dimension 
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iV = qK + 1 with K > 1. We feel that partially nonlocal models with nonva- 
nishing elementary lengths might find a very natural area of applicability in 
quantum graphs since the experimental waveguides and other nanotechnolog- 
ical realizations of quantum graphs almost certainly contain an uncertainty 
in the localization related to the degree of idealization of the real physical 
system in question [48] . In addition, the observability of the coordinate in 
a quantum graph may prove overridden by the transfer of emphasis to some 
other measured quantities (cf., e.g., the famous question "Can one hear the 
shape of a graph?" as asked in Ref. [19]). 



5.1 Numerical proofs 

We believe that even the oversimplified discrete quantum graphs with not 
too large q and/or iV can offer a new source of entirely abstract elementary 
models with, say, an unusual or anomalous parameter-dependence of their 
spectra tractable by numerical techniques. Having this purely descriptive 
ambition in mind let us now study the first few q = 3 models in some detail, 
emphasizing that a key to all of the above-sampled applications of non-Dirac 
metrics 9^7 lies in the demonstration of the reality of the spectrum of the 
initial Hamiltonian H which is non-Hermitian, H ^ in H.( F \ 

5.1.1 The Y-shaped discrete quantum graph with N = 7 




Figure 1: The spectrum of H^(X). 
The first nontrivial discrete q = 3 graph 



£5 



0x 2 X X 3 <> 



x 6 



X\ 



(34) 



X4 
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leads to the seven-dimensional one-parametric family of Hamiltonians 
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(J) 
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-1 + A 
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(35) 



They exhibit a particularly tight mutual interaction between the endpoints. 
The analysis of the energy spectrum pertaining to H^ 7 \\) may rely on 
the construction of the secular polynomial which appears factorized into its 
quadratic and quintic component. Thus, two of the levels are prescribed by 
explicit formulae, i?2,5 = 2 =)= yl — A 2 , while the remaining ones follow from 
the reduced secular equation, 

E 5 - 11 E A + (A 2 + 43) E 3 - (7 A 2 + 72) E 2 + (14 A 2 + 48) E - 7 A 2 - 9 = . 

This confirms that the energy levels of our seven-point Y-shaped quantum 
graph remain real in the interval of couplings A G (—1, 1). Its endpoints 
coincide with the position of the Kato's "exceptional points", i.e., of the 
values at which the first merger and complexification of a pair of energies 
takes place. 

The overall A— dependence of energies is displayed in Figure [TJ We see 
there that the spectrum has four fragile (i.e., asymptotically complex) and 
three robust (i.e., never complexifying) components. This observation fits the 
pattern predicted by the generic tunable VT— symmetric model of Ref. 



5.1.2 The next, q = 3 model with N = 10 
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Figure 2: The spectrum of iJ^ 10 ^(A). 
On the ten-point graph-lattice 

~X8]<>X5 X 2 II Xq II X 3 Xq<)[x^ 



X4 



(36) 



x 7 
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our Hamiltonian H^ 10 \X) acquires the matrix form 
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The set of its eigenvalues comprises the constant and doubly degenerate 
doublet E 5i q = 2, the two explicit roots E 3t $ = 2 ± a/2 — A 2 and the six 
implicit nodal zeros of the reduced secular polynomial 

E 6 — 13E 5 + (A 2 + 63) E 4 - (9 A 2 + 140) E 3 + (25 A 2 + 14l) E 2 — 

- (22 A 2 + 56) £ + 5A 2 + 6 = 0. 

The A— dependence of these energies is displayed in Figure |2] where the thick- 
ness of the middle straight line emphasizes that the exceptional constant- 
energy level E = 2 is doubly degenerate. 



5.2 Nonnumerical proof 

The most straightforward rigorous proof of the reality of the energies [i.e., 
of the reality of the spectrum of Hamiltonian H^ N \X)) may proceed via 
the explicit constructive demonstration of existence of at least one metric 
= 0(H) ^ / which makes this Hamiltonian Hermitian in 7i^ s \ 



5.2.1 The local versions of the discrete quantum graphs 

For our Y-shaped graphs the dimension N = 3K + 1 is finite so that we 
may search for special solution Bo of Eq. (|27|) using a diagonal matrix ansatz 
and some computer-assisted symbolic manipulations. In this way we verified 
that at N = 7 the diagonal solution is positive definite and, up to an overall 
factor, unique, 



e (7) 

(diagonal) 



' 1 





" 







1 
1 
1 

















1 
J 



At any iV = 3K+1 with K — 3, 4, ... we then revealed that the verification of 
the absence of any non-diagonal elements in the difference G 



(7) 

(diagonal) 



I can 



be performed non-numerically. Finally, using the assumption of diagonality 
we reduced the matrix difference H^Q — O if in Eq. (|27l) to the mere pair 
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of equations which specified the last two missing matrix elements in our 
ultimate solution compatible with Eq. (|27[) at any integer K, 



e (3JC+i) 

(diagonal) 



1 
1 
















••• 



1-A 
1+A 

u 1-A 



(37) 



Obviously, this matrix is invertible, Hermitian and positive definite so that 
it may play the role of the metric inside the whole interval of couplings 
A G (—1, 1). This confirms that our Hamiltonian H^(X) becomes Hermitian 
in the ad hoc Hilbert space TC^ where the diagonal metric (I3T1) is employed. 
Thus, we may modify our notation, write 



(3K+1) 
(diagonal) 



0^ ) (A)andH( 5 ) = 

Ti.Q N \\) and re-read the latter statement as the rigorous proof of the reality 
of the energies for A G (—1, 1). 



5.2.2 Equivalent Hermitian Hamiltonians 

Our constructive proof of existence of the (unique) diagonal metric O = 
Oq (A) given by Eq. ( 1371) implies the survival of the observability of the RK 
coordinates in both our (unitarily equivalent) physical Hilbert spaces 
and, in an amended notation, = Hq N \X). As long as the diagonality 

of ©[/^(A) is specified in coordinate representation, the usual multiplicative 
operator of coordinates remains Hermitian in the same Hilbert space H,q N \X) 
as the Hamiltonian H^ N '(X), indeed. 

In such an exceptional case it makes sense to recollect the Dyson-mapping- 
related factorization O = Q^Q of our diagonal metric and to restrict our 
attention to the positive definite and diagonal operator factors f2 = v^O. 
Their knowledge enables us to recall definition ( 1251) and, for illustration, to 
evaluate the related Hermitian isospectral partner Hamiltonian f) = (^^(A), 
say, at iV = 7, 
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The block-tridiagonal generalization of this formula to all dimensions N = 
3K + 1 is obvious. 



6 Manifestly nonlocal quantum graphs 

Our general quantization recipe described in paragraph 14.21 admits the tran- 
sition from the diagonal metric Oq [exemplified by Eq. ( J37l) at q — 3] to its 
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arbitrary non-equivalent alternative (1281) . This means that in the spirit of 
trivial examples studied in Ref. |TJ we are allowed to violate the locality also 
in all the other VT— symmetric quantum graphs. Moreover, we can demand 
that the sequence of the not necessarily unique nondiagonal metrics 0i, G>2 
... is partially ordered with respect to their increasing degree of non-locality 
defined, in a way suggested in Ref. [I], as a suitable growing function 6 = 9j 
of subscript j which is proportional, say, to the number of nonzero diagonals 
in the matrices or metrics Qj. 

Expansion (1281) of each individual Qj combines, in principle, several in- 
definite pseudometrics V = V(H). In this sense, our main task is twofold: 
Firstly we have to find at least one solution of Eq. (129]) . the nonlocality of 
which saturates the number 9j. Secondly we must guarantee the positivity 
of the resulting multiparametric sum Qj. 



6.1 Sparse- matrix pseudometrics at N = 4 



For our model (1311) the brute-force solution of Eq. ( 1261) leads to the following 
most general and exhaustive five-parametric formula for the (pseudo)metric, 



eV>(a,b,f,j,k) 



b b-f-j 
b f 
b J 



a f 

b — f — k + a 
k 



J 
k 

b — j — k + a 



(3? 



The variability of one of the parameters is spurious since it merely signals 
the double degeneracy of one of the eigenvalues. The fifth degree of freedom 
may immediately be interpreted, therefore, as an inessential angle of rota- 
tion in the two-dimensional subspace spanned by the corresponding pair of 
eigenvectors. The remaining four real parameters are independent and their 
presence reflects the well known ambiguity of the assignment of the metric 
to a given Hamiltonian (for a thorough discussion of this mathematical 
subtlety with serious physical consequences cf., e.g., Ref. [T7]). 



6.1.1 Positive-definite cases (metrics) 

The apparent simplicity of formula (1381 is slightly misleading because the 
interpretation of the matrix Q^>(a,b, f,j,k) as a metric requires that we 
guarantee its positive definite status [T7]. At N = 4, this property would be 
equivalent to the positivity of all of its four eigenvalues Tj > 0, j = 0,1,2,3, 

r 01 = a + b/2±bVl3/2, 

r 2>3 = a + b-f-j-k± v 7 2 • J 2 ■ /■••' - /./ • //■• • /■-./ • (39) 

We see that the positivity of the metric (i.e., of the norm) is guaranteed by 
the specification of the allowed domain T> of our quintuplet of parameters. 
This is particularly easy when we restrict our attention to the subdomain of 
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T> where / = j = k = and where 



a b b b 

b b + a 

b b+a 

b b+a _ 

We obtain the complete positivity constraint 

2a > \b\ (v / 13 + 1) , 6<0, 

2a>&(Vl3-l), &>0. 

This means that the allowed values of b belong to an interval which grows 
with a > 0. As long as this guarantees the positivity of r i at all f,j, k, the 
specification of the allowed domain of parameters will be completed by the 
inequality 

f + j + k+ \ / f 2 +j 2 + k 2 -fj-fk-kj<a + b. (40) 

Out of the doublet of constraints r 2> 3 > this is equivalent to the stronger 
one. We see that neither of the three parameters /, j, k will be allowed to get 
too large in comparison with a. 



(4) (a,&, 0,0,0) = 



6.1.2 Indefinite cases (generalized parities V) 



Elementary N = 4 example looks particularly well suited for illustrative pur- 
poses. Thus, in a search for the simplest possible parity-type pseudometrics 
V^(a,b, f,j,k) we have to construct such a solution of Eq. ( 129]) which is 
invertible but which is not positive definite. This matrix may be represented 
by the same formula as the metric B*- 4 -* (a, b, /, j, k) but at least one of the pos- 
itivity constraints (|39|) must be violated. In such a setting the requirement of 
maximal simplicity may start from the elimination of b which is exceptional 
in occurring nine times in Eq. (1381) . At b = we may also normalize a = 1 
[P^(a, 0, f,j, k) wouldn't be invertible at a = 0] and have 



(1,0 J,j,k) 
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(41) 



For a maximal simplicity of this matrix we leave just one of its parameters 
nonzero and get, say, 



P (4) (1,0,0,0,1) 



" 1 








" 





1 

















1 








1 






(42) 



In VT— symmetric models this matrix can play the role of parity V. Geo- 
metrically, it realizes the left-right reflection of our Y-shaped graph (!30|) . 
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6.2 Block-tridiagonal pseudometrics 



6.2.1 Solutions of Eq. ((291) at iV = 7 and N = 10 



The main source of insight in the structure of metrics and pseudometrics lies 
in the natural partitioning of Hamiltonians in q— dimensional subma- 

trices. At q = 3 the first two nontrivial though still sufficiently sparse matrix 
solutions were obtained by the straightforward computer-assisted symbolic 
manipulations with Eq. fl29|) . The results 
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(43) 



and 
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(44) 



open the way toward extrapolations. 



6.2.2 Extrapolation to any iV = 3K + 1 

The knowledge of the nontrivial solutions (j4"3]) and (jUJ) of Eq. (I29I) inspires 
the proposal of the following block-partitioned ansatz 
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Using the Hamilotnian of Eq. (1211) in its q = 3 version 
3 I -1 -1 -1 I 
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■1 — A 







-1 + A 








2 






2 



(were only the right low corner is coupling-dependent) and performing the 
appropriate insertions one readily verifies that Eq. (|29|) becomes an identity 
provided only that the unknown submatrix d(X) is defined by the elementary 
formula d(X) = — c(A) which remains the same at all integers K = 2, 3, . . .. 

Naturally, our block-tridiagonal ansatz (14"5l) as well as its verification and 
subsequent conclusions may immediately be extended to the other star graphs 
with q = 4, 5, . . .. The details are left to the reader. In what follows we 
shall address, instead, the other two questions, viz., a transition from the 
block-tridiagonal pseudometrics (j4"5]) to their block-pentadiagonal and higher 
descendants (cf. paragraph 16.31 below) and a transition from the indefinite 
pseudometric matrices [exemplified here by Eq. ( 1451) ] to the acceptable and 
positive definite band-matrix metrics expressed by the first nonlocal formula 



= 



(AT) 



(N) 



(46) 



[cf. paragraph 16.41 below and note that the latter expression is just the two- 
term truncated version of the general expansion (T28]) ]. 



6.3 Block-pentadiagonal pseudometrics 

The appeal of finding a block-pentadiagonal pseudometric (denoted by the 
symbol here) would lie in its possible insertion in the next truncated 
version of formula (j28|) . 

Q = 02 ]= p0W +1 Vr+vi N K (47) 

This formula may be used to define the more smeared, block-pentadiagonal 
nonlocal metrics. 

Once we leave the positivity questions aside and choose N = 10, the 
application of the computer-assisted direct-solution algorithm produces the 
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pseudometric solution V\ (A) of Eq. (129]) in the form 
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In a way recommended at q = 2 in Ref. pQ) this q = 3 solution has been 
made unique by the requirement of having a minimum of nonvanishing ma- 
trix elements in the first row. In our present case the optimality of such 
a requirement is less obvious. Indeed, as long as we have to optimize the 
sum ( H7|) rather than its individual components we may feel dissatisfied by 
the comparatively high number (= 39) of nonvanishing matrix elements in 
7^ 10) (A) [and, among them, 12 manifestly A— dependent items]. In such a 
case we may contemplate V^i^) given by the formula 
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containing just a minimum - 30 pieces - of the nonvanishing matrix ele- 
ments. We may also ask for the absence of fractions at a cost of having 32 
nonvanishing matrix elements in V^°\\) = 
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Another option could be based of the compromising choice of V^r'iX) = 
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with some elementary fractions but with minimum (= 8 pieces) of A— dependent 
matrix elements. 



6.4 Metrics as positive definite superpositions of pseu- 
dometrics 

6.4.1 Positivity constraint at N = 7 

The results of preceding paragraph have to be complemented by the empirical 
observation that the candidate for the metrics which is chosen in the one- 
parametric block-tridiagonal form ©fgf = /?xGq +v[ N ^ need not necessarily 
be positive definite. This may numerically be confirmed not only at vanishing 
(5 — but also at the positive values of (3 which are not sufficiently large. 
For illustration we selected (3 = 1/10 and found that in dependence on the 
value of A, three or four eigenvalues of ©n/ioiiA) remained negative. 




X 



Figure 3: The three lowest eigenvalues of the matrix 2 G + V\ . 

(7) 

One must be careful even if the candidate matrix 0L (A) looks domi- 
nated by its diagonal and safely positive-definite metric component. This 
is illustrated in Figure [3] where we displayed the three lowest eigenvalues of 
0^ at /3 = 2. In the picture we also see that there already exists just single 
eigenvalue which breaks the positivity and stays negative inside the whole 
interval of A G (—1, 1). 
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(7) 

Figure 4: The spectrum of the metric Oj 3 j 



3 ei 7) + v\ 



.(7) 



Our subsequent Figure H] illustrates the situation in which (3 = 3 is suf- 
ficiently large. Similar pictures can offer a comparatively reliable graphical 
confirmation of the positivity of any candidate ( l28i) for the metric. Thus, in 
our particular illustration we see that for 0f5(A) considered in the whole in- 
terval Of A e (-A< numer<ca, >(/?), \(numerical)^ j t j g su ffi cient to choose (3 = 3. 

Then our picture also leads to the graphical estimate of \( numerical ) (3) ~ 1. 
As long as the dimension N = 7 is small, this estimate may be replaced by 
the rigorous identification of \( numerical ) (3) = 1. By means of elementary al- 
gebra it is easy to show that this value coincides not only with the singularity 
(i.e., with the point of divergence) of the maximal eigenvalue of 0fJ, (A) but 
also with the zero of the minimal eigenvalue of the same matrix. 




(7) 

Figure 5: The lowest two numerical eigenvalues of metric 0[ 5/ / 2 ](^) 

(7) 

The A— dependence of the minimal eigenvalue of 0L- (A) is of particular 
relevance in the regime where the value of /3 decreases below it value used 
in Figure HI In our next Figure [5] we use (3 = 5/2 and see that the loss of 

(7\ 

the positive-definiteness of 0Lj(A) may be expected to occur at A = where 
the lowest eigenvalue would vanish. Of course, as long as the dimension of 
our illustrative model is small, it is very easy to find the corresonding critical 
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value of 
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minimal 
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_ + _ V44 + 36iV107+ — 



3 y/te + 36 iVW ' 



(4£ 



This quantity lies, in rational arithmetics, inside interval (39/16,5/2) and is 
numerically approximated by ~ 2.46050487. In a few complementary tests 
we found that the positivity of the metrics @^j(A) is still reliably confirmed 
at (3 = 149/60 ~ 2.483333333 since in the standard precision of computer 
arithmetics the related minimum ~ 0.02 of the lowest numerical eigenvalue 
is still safely positive at A = 0. 




Figure 6: The lowest two numerical eigenvalues of matrix @| 3 g/ 16 ](A) 

(7) 

Matrix O^j(A) retains its applicability as a metric also for /5s which lie 
slightly below their universal, A— independent bound fj48l) . In these cases 
one must restrict the admissible variability of the coupling A to intervals 
A G (-1,-A( ad hoc \(3)) and A G (\( ad hoc \f3) , l) . For explicit numerical 
illustration of such a conditional, A— dependent positivity of the metric below 
the critical boundary (148]) we choose (3 = 39/16 = 2.4375 and revealed 
that the lowest numerical eigenvalue of B (with the minimum ~ —0.02 < 
at A = 0) remained negative in the interval of A G (-A (ad hoc \ A (ad /loc) ) 
where \^ ad hoc ^ ~ 0.5 at our sample value of (3 = 39/16. It is necessary to 
keep in mind that A^ hoc \(3) quickly converges to one with the decrease of 
P < ^minimal- This is well illustrated by our last Figure [6] where we obtained 
X (ad hoc) _ 94 &t p = 38 / 16 = 2.375. 
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6.4.2 The positivity constraint at iV = 3K + 1 > 10 

For the ten-dimensional metric candidates for the metric with the same block- 
tridiagonal structure, e.g., for (3 = 3 in 
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very similar results were obtained supporting the applicability of our above- 
presented considerations to all the sequence of metric candidates 0^(A) 
with sufficiently large parameters (3 > (3 m inimai(N) and with unconstrained 
dimensions N = 3K + 1 = 13, 16, . . .. 



7 Discussion and summary 

In a way emphasized by several authors [33, [51] one should, strictly speak- 
ing, distinguish between the x— dependence in the wave function ip(x) and 
the x— dependence in the potential V(x) since in these two functions the con- 
cept of locality has a different mathematical background as well as physical 
meaning. Usually, the variable x entering wave functions ip(x) is treated as 
a measurable (i.e., real) quantity while the choice of the local V{x) may be 
treated just as a very special case of its possible generalized, equally admis- 
sible non-local alternatives. 

The more widespread use of the "non-local" wave functions ip(x) (where x 
need not be an observable real coordinate) only occurred during the growth 
of popularity of differential-operator Hamiltonians H^ PT ) = p 2 + V(pt)(x) 
where x has been considered complex |41j. A not too dissimilar non- locality 
also characterizes our present nonlocal versions of VT — symmetric discrete 
quantum graphs where we left the physical meaning of the spatial coordi- 
nate unspecified, citing only the related thorough discussion of this question 
available in our preceding paper Ref. [I]. 

We may summarize our present results by saying that we transferred 
the concept of VT— symmetry to the class of quantum systems living on 
graphs. These graphs generalize the usual real line of coordinates in one 
dimension. Several non-Hermitian, VT — symmetric versions of these struc- 
tures have been studied. On technical level we found one of the most vital 
mathematical sources of encouragement in a few older papers [52] whose au- 
thors demonstrated the practical viability of an approximative reduction of 
the graph edges to discrete lattices of points and vice versa. On this back- 
ground we succeeded in combining the existing quantum-graph concepts with 
the very fresh formalism using VT— symmetric Hamiltonians which are only 
made Hermitian via a comparatively complicated ad hoc inner product. 



27 



The technical feasibility of such a synthesis had several independent rea- 
sons. First of all, the spectra of energies proved real for the range of couplings 
A which stayed independent of the changes of the dimension N of the lattice. 
Second, our choice of the model proved lucky in the sense that in the current 
coordinate basis one of the constructed metric matrices ^ I happened to 
remain strictly diagonal. Together with the elementary form of matrix ele- 
ments of this particular metric O = Oo this not quite expected result made 
the necessary rigorous proof of the reality of the energies virtually trivial. 

As the first byproduct of this circumstance one of the eligible physical 
interpretations of our apparently non-Hermitian quantum-graph system re- 
mained trivial in the sense that it just required an inessential modification 
of the concept of observables and that it enabled us to construct its spec- 
trally equivalent representation characterized by the Hamiltonian which is 
Hermitian in the current sense (cf. operator {) in paragraph 15.2.21 above). 

In the same theoretical framework the second important consequence of 
the existence of the well-defined interval of admissible couplings may be seen 
in the emergence of new freedom in the choice of alternative, different physi- 
cal interpretations of the same quantum-graph Hamiltonian H(X). We were, 
once more, lucky in revealing that there exists an extremely natural partial 
ordering of these interpretations dictated merely by the degree of their non- 
locality or, in other words, by the extent of the smearing of the coordinate 
(the degree of this smearing or, if you wish, fundamental length 9j) grew 
with the subscript j of the selected closed-form metric Qj). 

On descriptive side let us re-emphasize the minimality of our interactions 
which were not supported by the whole graph but just by the closest vicinity 
of its endpoints. This also contributed to the feasibility of our constructions 
for which we had to develop several computer-assisted auxiliary symbolic- 
manipulation techniques and adaptive algorithms. Fortunately, the explicit 
calculations which were performed at the smallest dimensions usually gen- 
erated the output which admitted an extrapolation. Hence, the subsequent 
adaptation of the algorithms often degenerated to the mere verification of 
the extrapolated ansatz. 

During these constructions we completely avoided the unnecessarily com- 
plicated direct construction of the non-diagonal Dyson-map matrices Q and 
restricted our attention just to the metrics. Moreover we revealed that these 
metrics can be decomposed into sums of certain sparse pseudometrics, i.e., 
matrices with a sufficiently large portion of matrix elements equal to zero. 
This facilitated our calculations at higher dimensions. 

Our requirement of a fixed nonlocality does not make the resulting met- 
ric O = 0(H) unique. The constructive analysis of this metric-ambiguity 
problem in the specific quantum-graph setting can be perceived as an exten- 
sion of several recent non-graph (or trivial-graph, q = 2) studies assigning 
several non-equivalent probabilistic interpretations to a given Hamiltonian 
|53j . A partial correspondence can be then seen to the standard transitions 
between the coordinate and momentum representations of wave functions 
ip(x) where the role of the (unitary) Fourier transformation of Hilbert space 
is being taken over by the manifestly non-unitary Dyson mapping Q. In such 
a setting we made use of the fact that the argument x of wave functions need 
not necessarily carry the direct physical meaning of an (arbitrarily precisely 
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observable) coordinate. In terms of measurements the immediate connection 
between the coordinate x and its observability is, therefore, weakened. The 
coordinates can be interpreted as "smeared" [53]. In the context of topologi- 
cally nontrivial quantum graphs the prospective utilization of such a feature 
of phenomenological models looks particularly promising. 

In the language of mathematics our present family of discrete quantum- 
graph models proved exceptionally friendly. Their choice enabled us to disen- 
tangle the hidden Hermiticity constraints (1261) and to find closed formulae for 
the sparse-matrix metrics. The resulting availability of their generic multi- 
parametric forms has been interpreted as a new freedom of a phenomenology- 
friendly choice among alternative inner products specifying the non-equivalent 
physical Hilbert spaces of states Hf^K J ' — 0, 1, — Whenever asked for, an 
extension of our present particular quantitative and illustrative results on 
VT— symmetric quantum graphs to the higher degrees of nonlocality and/or 
beyond their equilateral q— point-star special class with small q = 3,4,... 
looks comparatively easy and straightforward. 

In the context of physics our present results are unexpectedly encourag- 
ing. A new flexibility of the model-building has been achieved here, first of 
all, via extension of the class of eligible interactions and, secondly, via the 
related innovative control of a degree of nonlocality reflected by the intro- 
duction of the "tunable" inner products. A deeper investigation of these pos- 
sibilities seems to form a new and promising quasi-Hermitian-graph project 
filling a certain gap in the broader context of existing directions of the study 
of quantum theory on graphs. 
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